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I. INTRODUCTION 

A study of the (static) inner fission barrier heights Bj* 
of even-even nuclei is motivated by the importance of this 
quantity for several physical phenomena. Many heavy 
nuclei decay by spontaneous fission, and the size of the 
fission barrier is a measure for the stability of a nucleus 
reflected in the spontaneous fission lifetimes of these nu- 
clei [l| . The probability for the formation of a superheavy 
nucleus in a heavy-ion-fusion reaction is also directly con- 
nected to the height of its fission barrier Q . The height 
_B^* is a decisive quantity in the competition between 
neutron evaporation and fission of a compound nucleus 
in the process of its cooling. The large sensitivity of the 
cross section a for the synthesis of the fissioning nuclei 
on the barrier height i?^* stresses a need for accurate cal- 
culations of this value. For example, a change of Bj* by 
1 MeV changes the calculated survival probability of a 
synthesized nucleus by about one order of magnitude or 
even more • The population and survival of hyperde- 
formed states at high spin also depends on the fission 
barriers [3, 4]. In addition, the r— process of stellar nu- 
cleosynthesis depends (among other quantities such as 
masses and /3-decay rates) on the fission barriers of very 
neutron-rich nuclei [1, 0] . 

During the last decade the role of triaxiality in the 
region of the saddle point of fission barriers has been rec- 
ognized and tested in many theoretical frameworks. It 
was found that the height of the barrier is reduced when 
triaxial shapes are allowed [7, 8]. However, this lowering 
strongly depends on the proton and neutron numbers 
and on the model employed. The investigations of in- 
ner fission barriers with triaxiality included are available 
within the frameworks of the microscopic-|-macroscopic 
method I9l4l3|[. the extended Thomas-Fermi plus Strutin- 
sky integral |14| , and non-relativistic energy density func- 
tionals based on Skyrme QIIMl and Gogny [i,[il[ll] 
forces. 

Covariant density functional theory (CDFT) [20] is 
an approach alternative to the above mentioned non- 
relativistic methods. Built on Lorentz covariance and 



the Dirac equation, it provides a natural incorporation of 
spin degrees of freedom [2ll, |23 and an accurate descrip- 
tion of spin-orbit splittings [22| (see also Fig. 2 in Ref. 
[23}). which has an essential influence on the underlying 
shell structure. Lorentz covariance of the CDFT equa- 
tions leads to the fact that time-odd mean fields of this 
theory are determined as spatial components of Lorentz 
vectors and therefore coupled with the same constants as 
the time-like components [l^l which are fitted to ground 
state properties of finite nuclei. In addition, pseudo-spin 
symmetry finds a natural explanation in the relativistic 
framework [25| . Over the years a large variety of nuclear 
phenomena have been successfully described within the 
CDFT (see Ref. ^2d\ and references therein). 

However, the progress in the study of the fission bar- 
riers within CDFT has been slower than in its non- 
relativistic counterparts. Inner fission barriers in sev- 
eral nuclei have been calculated in the axially symmetric 
relativistic mean field (RMF) + BCS approach in Refs. 
[26l - [30| . However, these investigations employ the con- 
stant gap approximation in the BCS part. Our recent 
study of pairing schemes used for the calculations of fis- 
sion barriers clearly shows that this approximation leads 
to unphysical results for the fission barriers [3l| . Thus, 
the results of these works have to be treated with a cau- 
tion. Fission barriers have also been studied in axially 
symmetric RMF calculations within the BCS approxi- 
mation using an effective density-dependent zero-range 
force in the pairing channel; this force represents a much 
more realistic approximation for pairing 31]. Recently 
also relativistic Hartree-Bogoliubov (RHB) calculations 
with the Gogny force DIS and with 5-forces in the pairing 
channel have been carried out [Slj for a study of fission 
barriers with axial symmetry. 

Unfortunately, axially symmetric calculations cannot 
be directly compared with experimental data since, as 
has been shown in non-relativistic calculations 0, Q , the 
lowering of fission barriers due to triaxiality is significant 
and can reach 3-4 MeV in some nuclei. At present, no 
systematic studies of the effects of triaxial degrees of free- 
dom on the height of inner fission barriers are available 
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in the covariant density functional theory; this degree of 
freedom has only been studied in specific nuclei such as 
264Hs [H and 240pu ^ within the RMF+BCS approach 
as well as ^''°Pu [SS*! within the RHB approach. Thus, the 
main goal of the current manuscript is to perform a sys- 
tematic investigation of the inner fission barriers within 
the triaxial RMF+BCS approach, and for the first time 
to confront these important experimental quantities with 
CDFT in a systematic way. 

The manuscript is organized as follows. The triaxial 
RMF+BCS theory and its details related to the calcula- 
tions of fission barriers are discussed in Sec. [TTl Sec. IIIII 
is devoted to the analysis of the effects of the truncation 
of the basis in the particle-hole channel of the model. 
Truncation effects in the pairing channel are considered 
in Sec. IIVI The results of the calculation of the fission 
barriers, the role of triaxiality and the comparison with 
experiment are discussed in Sec. |Vl Finally, in Sec. IVI[ 
we report on calculations with others relativistic param- 
eter sets based on density dependent coupling constants 
and in Sec. I VII] summarize the results of our work. 



II. THEORETICAL FRAMEWORK AND THE 
DETAILS OF NUMERICAL CALCULATIONS 

The starting point of Covariant Density Functional 
Theory (CDFT) is a standard Lagrangian density 

£ = -0 {^{id - g^uj - gppf - eA) - m - g^a) tp 

+ lid<y f - l^,.^'" + (1) 

which contains nucleons described by the Dirac spinors 
■0 with the mass m and several effective mesons char- 
acterized by the quantum numbers of spin, parity, and 
isospin. They create effective fields in a Dirac equation, 
which corresponds to the Kohn-Sham equation ISq in the 
non-relativistic case. 

The Lagrangian ([1]) contains as parameters the meson 
masses rria, rriuj, and nip and the coupling constants g^, 
g^, and gp. e is the charge of the protons and it vanishes 
for neutrons. This model has first been introduced by 
Walecka [H, [S^ ■ It has turned out that surface proper- 
ties of finite nuclei cannot be described properly by this 
model. Therefore, Boguta and Bodmer [38] introduced 
a density dependence via a non-linear meson coupling 
replacing the term im^cr^ in Eq. ([1]) by 

U{a) = Imla^ +^-g^a^ + \g,a\ (2) 

If not specified otherwise, the calculations are performed 
with the NL3* parameterization of the RMF Lagrangian 
[39j shown in Table H] Apart from the fixed values for the 
masses m = 939 MeV, = 782.6 MeV and nip = 763 
MeV it contains six phenomenological parameters 



TABLE L Parameters of the effective interaction NL3* in the 
RMF Lagrangian 



Parameters of NL3* 

m = 939 (MeV) 
m„ = 502.5742 (MeV) g„ = 10.0944 
m„ = 782.600 (MeV) = 12.8065 

mp = 763.000 (MeV) Qp = 4.5748 

g2 = -10.8093 (fm^^) 

go, = -30.1486 



9(7, gui, 9 pi 92, and 173, which have been recently adjusted 
to the experimental data in finite nuclei [39] eliminating 
a few deficiencies of the well known older parameter set 
NL3 0. 

In the current investigation, the triaxial RMF+BCS 
model is used [41|]. The RMF-equations are solved and 
at each step of the iteration the BCS occupation proba- 
bilities are determined. These quantities are used in 
the calculation of densities, energies and new fields for 
the next step of the iteration. We use monopole pair- 
ing force with the strength parameters Gt for neutrons 
(r = n) and protons (r = p); this method is based on 
the residual interaction of the seniority model 'S] . 

We start with a pairing strength parameters G and 
solve in each step of the iteration the gap equation ^42] 

fc>0 " 

with Ek = \J (sfc — A)^ + A^, where are the eigenval- 
ues of the Dirac equation and the chemical potential A 
is determined by the average particle number. Then the 
occupation probabilities 

and the gap parameters 

A^G^Ufc-yfe (5) 

/£>0 

are determined in a self-consistent way. The pairing en- 
ergy is defined as 

-Epair = -A^UfcWfe, (6) 

The sum over k in Eqs. ([3]), (O and ^ run over all states 
in the pai ring window < -Ecutoff- 

In Ref. |43| empirical pairing gap parameters 

= ^r' = |l7^ MeV (7) 

have been determined by the systematic fit to experi- 
mental data on neutron and proton gaps in the normal 
deformed minimum. 
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TABLE II: The G?, Ga , G? and G? parameters [in MeV] for 
diflerent parameterizations of the RMF Lagrangian and cutoff 
energy -Ecutoff = 120 MeV. 



Force 


G? 


G^ 


Gl 




NL3* 


9.1 


6.4 


8.1 


10.0 


DD-PCl 


9.2 


5.4 


8.0 


11.4 


DD-ME2 


9.2 


5.8 


8.1 


11.2 



These empirical gap parameters form the basis for the 
definition of the strength parameters Gt in the current 
manuscript. Two procedures have been used: a) In the 
analysis of different truncation (Sec. IIII[) and pairing 
schemes fSec. lIVp . the values Gn{Z, N) and Gp{Z, N) are 
defined for each nucleus with neutron and proton number 
N and Z under study from the requirement that, in the 
normal deformed minimum, the calculated pairing gaps 
coincide with the empirical values, b) In the systematic 
calculations of potential energy surfaces and fission bar- 
riers in actinides the same procedure is used first for all 
even-even nuclei in the Z = 90 - 100 and TV - Z = 42 - 66 
ranges resulting in a set of the strengths GrJZ, N) and 
Gp{Z,N). Then, the following expressions |44| 



N - Z 

A-Gn^G'l -G'i — MeV 



A 

N - Z 
-rv^nv'^^L MeV 

A 



(8) 
(9) 



are used in the calculations. The parameters G", GJ, 
G\ and G\ defined by the least square fit to the set 
of the Gn{Z,N) and Gp{Z,N). Their values depend on 
the parameter set of the Lagrangian and they are given 
in Table |lll In this way we have strength parameters for 
the effective pairing interaction depending in a smooth 
way on the neutron and proton numbers and, because of 
the changing level density, the gap parameters derived 
from those values show fluctuations as a function of the 
particle numbers. 

The calculations are performed imposing constraints 
on the axial and triaxial mass quadrupole moments. The 
method of quadratic constraints uses a variation of the 
function 



{H)+ J2 C'2MQ2,)-q2^) 
11=0,2 



(10) 



where (H) is the total energy, and (Q2p) denotes the 
expectation values of the mass quadrupole operators 

Q20 = 2z^-x^-y^ (11) 
Q22 = ~ y2 (j2) 



In these equations, 52^1 is the constrained value of the 
multipole moment, and G2^ the corresponding stiffness 
constants 42]. 

Correlations beyond mean field can influence the cal- 
culated values of flssion barrier height and the excitation 



energies of the superdeformed minima associated with 
the fission isomer |45j . The inclusion of rotational corre- 
lations can be performed by a symmetry restoration (an- 
gular momentum projection) and that of vibrations by 
a mixing of mean field states corresponding to different 
shapes by the method of generator coordinates (GCM). 
So far such an investigation has been performed only for 
240py within the generator coordinate method based on 
Skyrme DFT under the restriction to axially symmetric 
shapes 45]. It was found that compared to the ground 
state, angular momentum projection lowers the (axial) 
inner barrier by about 0.6 MeV and the fission isomer by 
about 1 MeV. In addition, it was found in Ref. [1^ that 
the schematic rotational correction based on the Belyaev 
moment of inertia ^465 frequently used in the literature 
gives a reduction of the fission barrier height which is 
appreciably larger than the one due to angular momen- 
tum projection. Based on these results, no rotational 
corrections are taken into account in our calculations. 
A similar approach has been used in the very successful 
calculations of Ref. 
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In the current investigation we do not consider the 
outer fission barriers. However, this restriction has its 
own merits. The inner barriers are generally better 
measured than the outer ones, and they are certainly 
more important for the r— process, since they determine 
thresholds. Furthermore, spontaneous fission lifetimes 
tend to be dominated by the inner barrier, even if oc- 
casionally an outer barrier can have a crucial effect if it 
is wide enough. The consideration of only inner fission 
barriers allows us to restrict our calculations to refiec- 
tion symmetric shapes, because the odd-multipole defor- 
mations (octupole, etc.) do not play a role in the in- 
ner fission barrier of the actinides and superheavy nuclei 
[13, 113, [13, 113, El- However, the analysis of the sym- 
metric fission pathway, which is the lowest in energy in 
heavy actinide nuclei such as some Fm isotopes (see Refs. 
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), is also possible in the current framework. 



III. TRUNCATION EFFECTS IN THE 
PARTICLE-HOLE CHANNEL 

The RMF-I-BCS equations are solved in the basis of 
an anisotropic three-dimensional harmonic oscillator in 
Cartesian coordinates characterized by the deformation 
parameters /3o and 70 and the oscillator frequency hwo = 
41/1-1/3 MeV (see Refs. [HI, El for details). The de- 
formation parameters of the oscillator basis /3q and 70 
are selected to be close to expected values /32 and 7 of 
constrained solution; this improves the convergence and 
minimizes the computational time. The truncation of the 
basis is performed in such a way that all states belonging 
to the shells up to Np fermionic shells and Nb bosonic 
shells are taken into account. The computational time 
increases considerably with the increase of Np but it is 
much less dependent on Nb- Thus, special attention has 
been paid to the selection of Np of the basis to be chosen 
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for a systematic study of fission barriers for the nuclei of 
interest, which provides at the same time a reasonable 
numerical accuracy in the predictions of the physical ob- 
servables. 

The selection of the truncation scheme was guided by 
the detailed analysis of the convergence performed in axi- 
ally symmetric RMF+BCS and RHB calculations of Ref. 
[3l| . In this reference, extensive tests of numerical con- 
vergence have been performed in the spherical, normal- 
deformed and superdeformed (/32 0.7 — 1.0) minima in 
the RMF calculations without pairing on the example of 
the nuclei ^^^V and ^°H20 with Z = 120 and N = 184. 
Contrary to the previous studies of the convergence in the 
RMF framework which were based on the comparison of 
the Np and Nf + 2 results, the "exact" solution (extend- 
ing the calculations up to Np = 36, Nb — 36) has been 
defined. Then it was shown that the binding energies 
and inner fission barriers for Np — 20 and Nb — 20 were 
described with an accuracy of approximately 200 keV 
and 100 keV, respectively, as compared with the exact 
solution. Therefore, the systematic calculations in the 
present manuscript have been carried out with Np — 20 
and Nb = 20. This selection of the basis is in line with 
our previous convergence tests in different mass regions 
(see, for example, Ref. [3]) which clearly show that at 
large deformations full convergence of the binding en- 
ergies is reached at larger values of Np than at lower 
deformations. In addition, they show that larger sizes of 
the basis (larger Np and Nb values) are needed for the 
nuclei with larger proton Z and neutron N numbers (see 
Rcfs. [11,141,150]). 

Of course, as long as the same number of fermionic 
shells Np and the same deformation of the basis /?o is 
used, calculations with the axial code should give iden- 
tical results to those obtained with the triaxial code at 
7 = 0°. In fact for nuclei under study we find agreement 
with an accuracy of approximately 50 keV throughout 
the deformation range of interest, which is caused by 
small differences in the mesh points of the Gaussian in- 
tegrations for the matrix elements. As a result for axi- 
ally symmetric shapes, the fission barrier heights which 
depend on the relative energies of the saddle point and 
normal deformed minimum differ by less than 50 keV in 
these two calculations. 

Extensive convergence tests in axially symmetric 
Skyrme-Hartree-Fock calculations also show that a simi- 
lar size of the basis is needed (see Sec. IIB in Ref. 48] for 
more details) . The comparison of these convergence tests 
suggests that there is no big difference in the convergence 
of total energies as a function of the size of basis in the 
relativistic and non-relativistic approaches. 

Because the fermionic basis contains large and small 
components of the Dirac spinor, the diagonalization of 
the Dirac equation is by a factor of approximately 8 more 
time consuming than the corresponding Schrodinger 
equation in the non-relativistic case. As a consequence, 
triaxial RMF-I-BCS calculations are more computation- 
ally demanding than the ones performed in the triaxial 



O OEq.(13) 

30MeV 
60MeV 
120 MeV 




10 



12 14 16 18 

Number of fermionic shells 



20 



22 



FIG. 1: (Color online) Neutron and proton pairing strengths 
as a function of number Np of fermionic shells for different 
pairing schemes. Pairing schemes are indicated either by cut- 
off energy iScutoft or by Eq. (fT3)) . 
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FIG. 2: (Color online) Neutron and proton pairing energies in 
the normal deformed minimum as a function of number Nf 
of fermionic shells for different pairing schemes. 



Skyrme EDF with BCS approximation of Ref. [?]. This 
is also a reason why we treat the pairing channel in the 
present triaxial RMF calculations in the BCS approxi- 
mation despite the fact that triaxial cranked Relativistic 
Hartree+Bogoliubov approach has been developed in the 
end of nineties [Sll - fssj and successfully applied to the de- 
scription of rotational structures in the pairing regime in 
different mass regions [51-55J. The RMF-fBCS calcula- 
tions are less time-consuming than the RHB calculations. 
In addition, as follows from our experience of the calcu- 
lations in axially deformed RMF-f-BCS and RHB codes 
[sH . the RMF-(-BCS calculations are more stable (espe- 
cially, in the saddle point region) than the RHB calcula- 
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FIG. 3: (Color online) The dependence of additional binding 
due to pairing on the number Nf of major fermionic shells 
for different pairing schemes. 



tions. 



IV. TRUNCATION EFFECTS IN THE PAIRING 
CHANNEL 

It is rather customary to analyze the dependence of 
total binding energies (or other physical observables) on 
the truncation of basis (see Sec. IIIip . However, we were 
not able to find any detailed investigation where the im- 
pact of the size of the oscillator basis on the parameters 
of pairing in the BCS framework has been discussed in 
detail. Thus, we studied the dependence of the strength 
of the pairing interaction on the number Np of fermionic 
shells under the condition that the proton and neutron 
pairing gaps in the normal deformed minimum are fixed 
for all values of Np and cut-off energies i?cutoff • 
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FIG. 4: (Color online) The dependence of fission barrier in 
■^^®Pu on the pairing scheme. Solid and dashed lines are used 
for axially symmetric and triaxial solutions, respectively. 
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FIG. 5: (Color online) The same as in Fig.|4]but for ^^''Cf. 



RMF-I-BCS calculations in the normal-deformed mini- 
mum have been performed with several values of the cut- 
off energy E'cutofi in Eq. ([SJ, namely i?cutoff = 30, 60, and 
120 MeV. In addition, the prescription of Ref. [5g| (indi- 
cated as "Eq. ([231)" in the figures) has been used. This 
prescription introduces smooth energy-dependent cut-off 
weights [13] 



fk 



1 



1 + exp[(efc - - AEr)/tJ.r] 



(13) 



for the evaluation of the local pair density. In this equa- 
tion, Efc are the eigenvalues of the Dirac equation and the 
chemical potentials At of the proton (r = p) or neutron 
(r = n) subsystems are determined by the particle num- 
bers Nr- The cut-off parameters AEr and /ir = AEr /lO 
are chosen self-adjusting to the actual level density in 
the vicinity of the Fermi energy. AE^- is fixed from the 
condition that the sum of the cut-off weights includes ap- 
proximately one additional shell of single-particle states 
above the Fermi surface 



feef2. 



(14) 



In Eq. (iMl) . fi,- denotes the single-particle space used in 
the calculations. 

Figures [U [5] and [3] summarize the results of this study 
for the normal-deformed ground state in ^"^^Pu. One can 
see in Fig. [T] that the strengths of the pairing interac- 
tion depends not only on the cut-off energy i^cutoff but 
also on the number Np of fermionic shells employed in 
the calculations. This dependence is very weak for the 
prescription of Ref. 156] because here the effective pair- 
ing window is quite small being around 7 MeV. On the 
other hand, the dependence of the pairing strength on 
Np increases with the increase of ^cutoff- This can be un- 
derstood in the following way: an increase of Np brings 
more single-particle states into the pairing window thus 
effectively requiring the decrease of pairing strength in 
order to keep the pairing gap fixed. This effect becomes 
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|32cos(7+30) 

FIG. 6: (Color online) Potential energy surface in ^*°Pu. The 
energy difference between two neighboring equipotential lines 
is equal to 0.5 MeV. The blue dashed line with solid circles 
shows the lowest in energy solution as a function of /32. Fur- 
ther details are given in the text. 



more pronounced for larger pairing windows, which ex- 
plains the steeper decrease of the pairing strength as a 
function of Np with increasing i^cutoff- 

The dependence of proton and neutron pairing ener- 
gies ii'paii. and -Epaij. on the cut-off energy E'cutoff and on 
the number Np of fermionic shells employed in the cal- 
culations is shown in Fig. [5] These energies depend only 
weakly on Np in the case of prescription of Ref. [s^ be- 
cause of the small effective pairing window. However, 
similar to the pairing strengths the dependence of pair- 
ing energies on Np increases significantly with increasing 
^'cutoff ■ The origin of this feature is the same as in the 
case of the pairing strengths; it is discussed above. 

Note, however, that the dramatic changes in the pair- 
ing energy cannot be seen directly in the change of the 
energy, because they are compensated to some extent 
by the fact that larger pairing seen in pairing energies 
causes a wider distribution of the occupation probabili- 
ties around the Fermi surface. Therefore, we study in 
Fig. [3] the energy difference between the binding energies 
^pair _ ^unpair obtained in two self-consistent calcula- 
tions with and without pairing. It turns out that for 
Np > 14 this difference, that reflects the real physical 
impact of pairing, is smaller than 1 MeV and it does nei- 
ther depend on the cut-off energy E'cutoff nor on the value 
of Np. Somewhat different values of EP''" - E'^'^p'^" at 
lower Np values are due to the fact that at these values of 
Np the effects of the truncation of basis in the particle- 
hole channel have not been eliminated (see Sect. IIIII for 
detail) . 

In Figs, m and [5] we compare the deformation energy 
curves for three different pairing schemes for the nuclei 



^^^Pu and ^^°Cf. The deformation energy curve for the 
axially symmetric solution is obtained as the 7 = 0° 
cross-section of the potential energy surface. The defor- 
mation energy curve for the triaxial solution is obtained 
by the minimization of potential energy surface along 
the /32-direction. We show the deformation energy curve 
for the triaxial solution only in the range of values 
where it is lower in energy than the deformation energy 
curve of the axially symmetric solution. Note that the 
potential energy surfaces are normalized to zero at the 
normal-deformed minimum. As discussed in Ref. [Slj we 
can see that these different schemes predict somewhat 
different fission barriers. In the systematic calculations 
presented in the following sections, we use the cut-off en- 
ergy -Ecutoff = 120 MeV. The selection of this value is 
based on the results of Ref. [3l| , where it was shown that 
the difference in the height of fission barriers obtained in 
the RHB calculations with finite range DIS Gogny force 
and zero-range (5-force is minimal when the cut-off energy 
EcntoS = 120 MeV is used (see Fig. 6 in Ref. [U). 



V. SYSTEMATIC ANALYSIS OF THE INNER 
FISSION BARRIERS 

In this section we carry out a systematic investigation 
of fission barriers of even-even nuclei in the actinide re- 
gion based on the parameter set NL3* and the pairing 
strength parameters given in Table HIl In Fig. [5] we show 
as an example the potential energy surface of the nu- 
cleus ^^"^Pu in the /3-7 plane. For axial symmetry we find 
the normal deformed minimum of the ground state at a 
deformation /? ~ 0.28, a maximum at /? ~ 0.52 and a 
superdeformed minimum at (3 ^ 0.96. We observe that 
the fission path (the part of blue dashed line between 
normal and superdeformed minima) bypasses the axial 
barrier between the normal and superdeformed minima. 
The barrier height is determined by the maximum of the 
energy along this fission path. 

The deformation energy curves for other even-even nu- 
clei in this region obtained in these calculations are shown 
in Fig. [T] Full black lines show axially symmetric solu- 
tions, while we show the values of the deformation energy 
curves along the triaxial fission path by red full curves. 
One can see that by allowing for triaxial deformation the 
fission barrier heights are reduced by 1 — 4 MeV as com- 
pared with axially symmetric solutions. This lowering 
depends on the proton and neutron numbers. It also 
brings in average the results of the calculations in closer 
agreement with experimental data shown by green solid 
circles in Fig.[7l These circles display the height of the ex- 
perimental fission barrier at the calculated /3-deformation 
of the saddle point. The calculated 7-deformations of the 
triaxial parts of the fission path are shown in Fig. [51 On 
average they are close to 10°. 

The microscopic origin of the lowering of the barrier 
due to triaxiality can be traced back to the changes of the 
level density in the vicinity of the Fermi level induced by 
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FIG. 7: (Color online) Deformation energy curves of even-even actinide nuclei obtained in the RMF+BCS calculations with 
the NL3* parameterization. Black solid lines display the deformation energy curves for the axially symmetric solution, while 
red solid lines show the deformation energy curves along the triaxial part of the fission path. Solid green circles show the 
experimental values of the height of inner fission barrier. Experimental data are taken from Table IV in Ref. [31 • A typical 
uncertainty in the experimental values, as suggested by the differences among various compilations, is of the order of ±0.5 
MeV [3. 
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FIG. 8: (Color online) The 7-deformations for the results of calculations shown by red lines in Fig. [T] 



triaxiality. Fig. |9] shows the Nilsson diagrams for protons 
and neutrons for the axially symmetric solution in ^^^Pu. 
The blue boxes in these diagrams define the deformation 
and energy ranges in which the axially symmetric and 
triaxial solutions are compared in Fig. 1101 The lower 
(upper) deformation in these boxes corresponds to the 
deformation range over which the triaxial solution (red 
curve in Fig. [7]) is lower in energy than the corresponding 
axially symmetric solution (black curve in Fig. [7]). The 
lower and upper energy values in these boxes are defined 
approximately as ep ± 3 MeV. 

Proton and neutron single-particle energies within 
these deformation and energy ranges are shown for ax- 
ially symmetric and triaxial solutions in Fig. 1101 One 
can see that the single-particle level density at the Fermi 
level is lower for triaxial solutions than for axially sym- 
metric solutions. This is especially clear at the deforma- 
tion corresponding to the saddle point of the axially sym- 
metric solution (indicated by vertical dotted blue lines 
in Fig. 1101) which correspond to a maximal level den- 
sity and maximal pairing correlations. A lowering of the 
level density at the Fermi surface leads to a more neg- 



ative shell correction energy (as compared with axially 
symmetric solution), and, as a consequence, to a lower 
fission barrier. This is in agreement with the analysis of 
Ref. which also attributes the lowering of the inner 
fission barrier due to triaxiality to microscopic (shell cor- 
rection) part of the macroscopic -l-microscopic model. A 
similar mechanism is responsible for the lowering of the 
asymmetric saddle with respect to symmetric saddle at 
outer fission barrier (see Sect. VI in Ref. [13 )■ 

Figures [TT] and [12] show the differences between calcu- 
lated and experimental heights of inner fission barriers. 
The average deviation between theory and experiment is 
0.76 MeV. This is comparable with the results obtained 
in the macroscopic-|-microscopic method (see Sec. IVC 
and Fig. 11 in Ref. [13 and Sec. VII A in Ref. [13) 
which describe experimental fission barriers with an av- 
erage error of around 1 MeV. 

It is necessary, however, to say that neither proton nor 
neutron particle number dependences of fission barrier 
height are completely reproduced in these calculations. 
This is clearly seen in Figs. [TT] and [T^ However, the same 
problem exists also in macroscopic+microscopic calcula- 
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FIG. 9: (Color online) Proton and neutron single-particle energies in ^*^Pu as a function of the quadrupole deformation P2 for 
the axially symmetric solution. Solid (black) and dashed (red) lines are used for positive and negative parity states, respectively. 
Blue boxes show the regions which are displayed in more details in Fig. llOl below. Fermi energies ep are shown by dot-dot-dashed 
(green) lines. 



tions (see Fig. 11 in Ref. [llj] and Figs. 23-32 in Ref. 
[1^). There are very few energy density functional cal- 
culations of the fission barriers with triaxiality included, 
and neither of them confronts in a systematic way experi- 
mental data in actinides. However, limited results in the 
Skyrme EDF presented in Ref. [16l | show similar unre- 
solved particle number dependences for the inner fission 
barrier heights. 



VI. RESULTS FOR THE PARAMETER SETS 
DD-ME2 AND DD-PCl 

In order to investigate to what extent our results de- 
pend on the density functional under investigation we 
performed also an analysis of the fission barriers of the 
two nuclei ^'"^Pu and ^36^ using the parameter sets DD- 
ME2 m and DD-PCl [59*1. The first is a representative 
of the class of the RMF models [58, 6(| where the nucleus 
is described as a system of Dirac nucleons interacting via 
the exchange of mesons with finite masses leading to in- 
teractions of finite range. An explicit density dependence 
for the meson-nucleon vertices is used. The DD-PCl pa- 
rameterization belongs to the class of the RMF models 



in which the finite-range meson exchange is replaced by 
zero-range interactions with density dependent coupling 
constants and derivative terms jsol . IsTI . [63 [ . 

In Fig. [T2] we compare the deformation energy curves 
of ^'"'Pu and ^se^j obtained in the calculations with the 
three parameter sets NL3*, DD-PCl and DD-ME2 of the 
RMF Lagrangian. Although there are some differences 
between the deformation energy curves obtained in the 
calculations with different parameterizations, in general, 
they show the same features. In addition, calculated fis- 
sion barrier heights reasonably agree with experimental 
data. More systematic investigations of fission barriers 
with the DD-ME2 and DD-PCl parameterizations of the 
RMF Lagrangian are in progress and their results will be 
presented in a forthcoming manuscript. 



VII. CONCLUSIONS 

We presented here the first systematic investigation of 
triaxial fission barriers in the actinide region within co- 
variant density functional theory. The calculations have 
been carried out with the parameter set NL3* and they 
have been compared in specific cases also with the results 
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FIG. 10: (Color online) A comparison of proton and neutron single-particle energies at axially symmetric (left panels) and 
triaxial (right panels) solutions. The same notation for the lines as in Fig. [9] is used. Vertical blue dotted lines are used to show 
the deformation of the saddle point obtained in the axially symmetric solution. 
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FIG. 11: (Color online) The difference between experimental 
and calculated heights of inner fission barriers as a function 
of neutron number A*'. 



of parameter sets DD-ME2 and DD-PCl. Pairing corre- 
lations are taken into account in the BCS approximation 
using seniority zero forces adjusted to empirical values of 
the gap parameters. It is found that with only one ex- 
ception (^•^■^Th) in all the nuclei under investigation the 
height of the inner fission barrier is reduced by allowing 
for triaxial deformations by 1 — 4 MeV. The fission path 
avoids a maximum of the axially symmetric potential en- 
ergy surface between the first and the second minimum 
by going through a valley in the (/?, 7) plane with a tri- 
axial deformation 7 w 10°. A systematic comparison of 
our results with experimentally determined fission barri- 
ers in this region shows reasonable agreement with data 
comparable with the best macroscopic -l-microscopic cal- 
culations. 
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FIG. 12: (Color online) The same as in Fig. [TT] but as a 
function of proton number Z. 
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